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SUPERCONVERGENCE POINTS OF FRACTIONAL SPECTRAL 

INTERPOLATION * 

XU AN ZHAO t AND ZHIMIN ZHANG* 


Abstract. We investigate superconvergence properties of the spectral interpolation involving 
fractional derivatives. Our interest in this super convergence problem is, in fact, twofold: when 
interpolating function values, we identify the points at which fractional derivatives of the interpolant 
super converge; when interpolating fractional derivatives, we locate those points where function values 
of the interpolant superconverge. For the former case, we apply various Legendre polynomials as basis 
functions and obtain the superconvergence points, which naturally unify the superconvergence points 
for the first order derivative presented in [Z. Zhang, SIAM J. Numer. Anal., 50 (2012), 2966-2985], 
depending on orders of fractional derivatives. While for the latter case, we utilize Petrov-Galerkin 
method based on generalized Jacobi functions (GJF) [S. Chen et al., arXiv: 1407. 8303vl] and locate 
the superconvergence points both for function values and fractional derivatives. Numerical examples 
are provided to verify the analysis of superconvergence points for each case. 

Key words. Superconvergence, fractional derivative, spectral collocation, Petrov-Galerkin, 
generalised Jacobi functions 

AMS subject classifications. 65N35, 65M15, 26A33, 41A05, 41A10 


1. Introduction. Superconvergence of numerical methods usually happens when 
the convergence rate at some special points is higher than the theoretical global 
rate [6ll9ll24M31j . The investigation of superconvergence provides a fundamental insight 
for post-processing and adaptive algorithm design, which leads to higher accuracy 
and more efficient numerical methods m- 

For classical problems, the superconvergence phenomenon of the h-version finite 
element method is well understood and addressed in the literature, see, e.g., mm, 
and references therein. On the other hand, there has been only a limited study of 
the superconvergence phenomenon for spectral methods [3TH33I . Specifically, in 
mm, Zhang identified the derivative superconvergence points of the Legendre and 
Chebyshev spectral collocation methods for the two-point boundary value problem. 
In addition, Zhang |33l investigated various polynomial interpolations and located 
superconvergence points for the function value and the first derivative. 

In this paper, we conduct the first superconvergence study of fractional operators. 
Since fractional calculus generalizes the classical (integer-order) differentiation and 
integration to any order, the study of superconvergent fractional derivatives unifies 
the theoretical investigation in l33l for the first-order derivative. 

Over the last two decades, fractional differential equations (FDEs) have been 
demonstrated to be more effective, when modeling some complex systems in physics, 
finance and et al. mmmmm- compared to classical models. These new models 
are usually derived by replacing integer-order derivatives with fractional derivatives 
in classical models. Meanwhile, there has been a growing need for the development of 
high-order numerical algorithms for solving FDEs. Due to the nonlocal definition of 
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fractional derivative, existing methods including finite difference and finite element 
methods [T7J[2ni[lll[2Sl[lil[17J[331l35] mostly lead to low-order schemes. 

Spectral methods are promising candidates for solving FDEs since their global 
nature fits well with the nonlocal definition of fractional operators. Using integer- 
order orthogonal polynomials as basis functions, spectral methods [TTlll2U26ll30| really 
help with the alleviation of the memory cost for discretization of fractional derivatives. 
Furthermore, to deal with singularities, which usually appear in fractional problems, 
in the works 00I2SI55] the authors design suitable bases. 

A major difficulty in the investigation of superconvergence of spectral methods 
for fractional problems, compared with integer-order derivatives, is the nonlocality of 
the fractional operator and the complicated form of fractional derivatives. The second 
challenge is the construction of a good basis for a spectral scheme. Given a suitable 
basis, one can then begin the analysis of the approximation error in order to locate 
the superconvergence points. 

In this paper, we interpolate function values using Legendre polynomials, and use 
the error equation to obtain superconvergent points for Riemann-Liouville fractional 
derivatives. It turns out that these superconvergence points are zeros of the frac¬ 
tional derivatives of the corresponding Legendre polynomials. We also consider the 
interpolation of fractional Riemann-Liouville derivatives, and apply a GJF-spectral- 
Petrov-Galerkin scheme [3] to solve an equivalent fractional initial-value problem. We 
have found that superconvergence points for numerical function values and fractional 
Riemann-Liouville derivatives are zeros of corresponding basis functions and Gauss 
points, respectively. 

The organization of this paper is as follows. In Section 2, we define notation and 
give properties of fractional derivatives and of GJFs. The interpolations of the func¬ 
tion values using various Legendre polynomials are presented in Section 3. Moreover, 
we analyse the error to identify superconvergence points for fractional derivatives of 
the interpolant. In Section 4, we show the locations of the superconvergence points 
both for the numerical solution and the fractional derivatives using a GJF-spectral- 
Petrov-Galerkin scheme. The numerical tests for both cases are displayed at the end 
of Section 3 and Section 4 respectively. We conclude with a discussion of our results 
and their applications in Section 5. 

2. Preliminaries. In this section, we present some notations and lemmas which 
will be used in the following sections. 

Definition 2.1. The fractional integral of order /i £ (0,1) for function f(x) is 


defined as 



( 2 . 1 ) 


Definition 2.2. Caputo fractional derivative of order p £ (0,1) for function 
/( x) is defined as 



Definition 2.3. Left Riemann-Liouville fractional derivative of order fi £ (0,1) 
for function f(x) is defined as 



x > x L , (2.3) 
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Definition 2.4. Right Riemann-Liouville fractional derivative of order fi G (0,1) 
for function f(x) is defined as 


( RL x VZJ)(x) 


1 

r(i-/x) 




f(s) 

(s — x)v 


ds, 


X < XL, 


(2.4) 


For the convenience, if p. < 0, we denote ( RR/ D£f)(x) = ( Xi I“ M /)(a;) throughout 
the paper. 

Lemma 2.5. (see w Let pgl and n-l<p<nel If f(t) and <p(t) along 
with all its derivatives are continuous in [a, t\. Under this condition the Leibniz rule 
for Caputo fractional differentiation takes the form: 

OO / \ 77* 1 . _ p 

c a vf mm] = e (l )v w (t) C a v t~ k f(t) - E 

k =0 ' ' k =o 1 ' 


Lemma 2.6. (see m Chap. 2) Let pel and n—l<p<n£R. If f(t) 
and ip(t) along with all its derivatives are continuous in [a,t\. Under this condition 
the Leibniz rule for Riemann-Liouville fractional differentiation takes the form: 

RL a v ? mm\ = e RL a vr k f(t). 


We recall the definition of GJFs, which extend the parameters of classical Jacobi 
polynomials to a wider range . 

Definition 2.7. (see ffff) [Generalized Jacobi functions] Define 

+ J$- a ’ p) (x) := (1 - x) a P^ p \x), fora > -1, /3 G R, (2.5) 

and 

~Jk a ’~ P) (x) := (1 +xfP[ a ^ ) {x), fora G R, /3 > -1, (2.6) 

for x G A = (— 1,1) and n Js 0. 

In what follows, we present the selected two special cases for the fractional deriva¬ 
tives of GJFs: 

• Let a > 0, /3 G R and n G No 


RL x V([ + j[ l -^\x)) 

_ r(n + a -(- 1) p(o, a + 0 )/\ 

n\ n 

(2.7) 

• Let /3 > 0, a G R and n G No 




_ T(n + 0 + l) p(a+0<o) ^ 

(2.8) 

For convenience , we introduce two sets of parameters for GJFs: 


+ TE :={(«,/3) :a>0, p> 

-1}, 

(2.9) 

+ T“’' 9 := {(a,/3) : a > 0, —a - 

- 1 < /3 = -k < -1, k G N}. 

(2.10) 


For (a, fi) G we define the finite-dimensional fractional-polynomial space: 

+ Er a,/3) ( A ) : = spanl+J^’^, 0 < n < N}. 


( 2 . 11 ) 
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3. Legendre interpolant. In this section, we interpolate a smooth function u 
at a set of TV + 1 special points , which are zeros of some orthogonal polynomials on 
[—1,1], In particular, we want to find un such that 


UN(x k )=u(Xk), -1 < X 0 < X\ < ■■■ < XN < 1, (3.1) 

where the interpolation points {xk, k = 0,1, • • • , TV} are zeros of say, the Legendre 
polynomials. Our goal is to identify those points y'p where ^V^un, (3 € (0,1), the 
fractional Riemann-Liouville derivative of the interpolant, superconverges to 
in the sense that 

TV a |“2?f(n-n n )(yf ) )| <C max \ R _ L {D^(u - u n ){x)l a > 0. (3.2) 


According to the description (13.21) of superconvergence points originated from [33], 
if (3 = 1 and {z/'-} are independent of the particular choice of u , we say that {z/'} are 
the superconvergence points for the first derivative of the interpolant. To step further, 
here we want to seek out the points y^ independent of the function u, at which the 
corresponding /3-th Riemann-Liouville derivative ^{D^un superconverges to ^V^u. 
In other words, for different order of fractional derivative, different superconvergence 
points are needed. 

3.1. Analysis. Considering polynomial interpolation for (13.11) . the key rule for 
locating the superconvergece points is to analyze the interpolation error. We are now 
in the position to show the theoretical analysis. 

Proposition 3.1. For the interpolation of (EED using collocation points as 
the zeros of Legendre polynomials Ln + \{x), Legendre-Lobatto polynomials (Lat+i — 
Lat_i)(s), and Legendre-Radau (right and left) polynomials (Ln + i ±Ln)(x), the /3-th 
Riemann-Liouville fractional derivative superconverges at {T^' 9 ' 1 }, which satisfy 

rl 1 VP Wn+ r(^ } ) = 0, k = 0,1,2, • ■ • , TV, (3.3) 


where wn+i(x) denotes the aforementioned four sets of polynomials respectively. 

Proof. Let u be analytic on / = [—1,1]. According to [33], u can be analytically 
extended to B p , which is enclosed by an ellipse E p with ±1 as foci and p > 1 as the 
sum of its semimajor and semiminor: 

E p : 0 = ^(pe i9 + p-'e- 16 ), 0 < 6 < 2tt. (3.4) 

We consider polynomial un G Pn who interpolates u at TV+1 points — 1 < Xo < X\ < 
■ ■ ■ < Xn < 1. The error equation is, according to [5], expressed as 


u{x ) — Un(x ) 


1 f w N+ i(x) u{z ) 
2ni J Ep wn+i(z) z - x 


(3.5) 


where wn+i(x) = c J}‘^L 0 (x — Xj). 

Noticing Lemma 12.61 taking the fractional derivative of the above equation we 
have 


-u n (x)] 


1 

27TZ 



' “ f(3\ ^- fc Wjv+1 Qr) 

few (Z-XY + 1 


u{z) 

w N +i(z) 


dz. 


(3.6) 
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Let us examine the error equation (13.61) . According to the analysis in [33], the 
leading term of the error is determined by ^{D^w n+\ (x). At the N special points 

= 0, • • • ,N}, which satisfy r 3)T>xWn+i{£,\ 0 ' > ) = 0: we have the remaining 
terms, which are usually smaller than the first term in magnitude at least by a factor 
N a . Here we interpolate the function using Legendre polynomials, Legendre-Lobatto 
polynomials, and Legendre-Radau (right and left) polynomials, thus the correspond¬ 
ing superconvergence points are zeros of RT {'D^ : vj n+1 (x). where wn + i(x) denote the 
aforementioned four sets of polynomials respectively. □ 

From Lemma [2.51 the superconvergence points for the Caputo fractional derivative 
is easily drawn. 

Remark 3.2. Similar to the proof in Propsition \3. il for the interpolation of 
(E3D using collocation points as the zeros of aforementioned four kinds of Legendre 
polynomials, the /3-th Caputo fractional derivative superconverges at {z^}, which 
satisfy 

SiD^WN+iizjP) = 0, k = 0,1,2, •. • , N, (3.7) 

where wn+\{x) denotes the aforementioned four sets of polynomials respectively. 

As an extension, superconvergence results are derived similarly for Chebyshev 
polynomials in the following. 

Proposition 3.3. For the interpolation of JUD using collocation points as 
the zeros of Chebyshev polynomials of the first kind Tn+i(x), Chebyshev-Lobatto 
polynomials (Tjv+i — Tn-i){x) , and Chebyshev-Radau (right and left) polynomials 
(T/v+i ± Tjv)(x), the fi-th Riemann-Liouville fractional derivative superconverges at 
{y^}, which satisfy 


“2^«i JV+1 (^ ) )=0,fc = 0,1,2, — ,N , (3.8) 

where wn+ i(x) denotes the aforementioned four sets of Chebyshev polynomials re¬ 
spectively. 

In what follows, we explain how to compute {£[, /3 '} hi (13.31) . 

Denote L n (x), the Legendre polynomial of degree n. In spectral methods, we often 
use the following four combinations. For convenience, express them as generalized 
Jacobi functions: 


L n (x) = Pi°’°\ x), (Legendre polynomial) 

(L n — L n - 2 ){x) = P( L ~ 1 ' 1 \x), (Lobatto polynomial) 
(L n + L n _ 1 )(x) = P^ (x), (Left Radau polynomial) 


(L n — L n -i)(x) = Pb 1,0 ' ) (x). (Right Radau polynomial) 

Using formula [10] 

p( a <fi( x ) - r ( n + <*-& + l)r(n + !) / 1 + aA 0 p( a -ff), s f 3 gi 

r n -p r(n _^ + 1)r(n + a + 1) ^ 2 J ^ W, ( 3 . 9 ) 
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and 


RL 

-l 


DtP: 




0*0 


r(n + a + fi + 1) (a+^,/ 3 +^ 1 ), , 
2^r(n + a + 1) 1 


we can derive 



2~^r(n + 1) 
T(n - /x + 1) 



(3.10) 


(3.11) 



(n + a)r(rt) ^ 1 + s V M ( a+A1 ,i_ M ), , 
2* i r(n — /x +1) V 2 / n_1 ’ 


Similarly, we have 


i?L 

X 




2-^r(n + l) 
T(n - n + 1) 



a = —1, 0. 

(3.12) 


(3.13) 


i?L 




(n + ct)r(n) 
2^r(n — /x + 1) 



i-/x 


p(l-/i,a+/i) 

-*n—1 


0 ), 


a = —1, 0. 

(3.14) 


Therefore, we conclude that m is equivalent to find zeros of: 

(1) -Pjv+i^ 0*0 f° r the Legenrde interpolation, i.e., interpolation at roots of Pjy+i 0*0 > 

(2) (1 + x) 1 ~ /J ’P^‘~ 1,1 ~ M \x) for the Lobatto interpolation, i.e., interpolation at 

roots of py+r 1 ^); 

(3) (1 +xY-^’ 1 ^(x) for the Left-Radau interpolation, i.e., interpolation at 
roots of P^y) (:r). 

Similar results can be obtained for the right fractional derivative RR D^. 

We see that a fractional derivative of a polynomial of degree N + 1 may have 
N + 1 roots (comparing with N roots for an integer derivative). 


3.2. Numerical validation. We report on the numerical results using the afore¬ 
mentioned four sets of Gauss interpolation points in (13.11) . 

Example 3.1 We choose u{x) = j^q(x + l) 1015 in (13.11) . 

Fractional derivative errors of its interpolant, using Gauss points, Gauss-Lobatto 
points, and the left Gauss-Radau points with number of points of IV = 12 are depicted 
in Figures 13.1113.31 respectively. From Proposition 13. II we know that for a polynomial 
interpolant of degree TV, the leading error term is a polynomial of degree N + 1, and 
there are N + 1 roots for its fractional derivative. From Figures 10113.31 it is clearly 
seen that there are 13 superconvergence points. Particularly, —1 is a superconvergece 
points for the Gauss-Lobatto and left Gauss-Radau interpolants. To view those roots 
near —1, we add zoom-in windows into these figures. 

Five different colors of asterisks denoting the superconvergence points, with orders 
/? = 0.1, 0.3, 0.5, 0.7, 0.9 for Legendre polynomial interpolation, are show in Figure 
ELD It is easily shown that for different orders of fractional derivatives, the errors 
at the superconvergence points are significantly smaller than the maximum error. In 
addition, it is also indicated that when the fractional order becomes smaller, the error 
curve goes lower. 
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Similar situations is easily drawn from Figure 1572113.31 for Legendre-Lobatto poly¬ 
nomial and the left Legendre-Radau polynomial respectively. While for approxima¬ 
tions using Gauss points, the errors on the boundaries are larger than that using 
Gauss-Lobatto points and left Gauss-Radau points. 



Fig. 3.1. Curves of ^_\T> x {u — un) for Example 3.1 using Gauss points, where * denotes the 
corresponding superconvergence points. 


4. Fractional derivative interpolation. In this section, we present the inter¬ 
polation for the fractional Riemann-Liouville derivatives of a smooth function using 
GJFs as basis functions and investigate superconvergence points for the function value 
approximation and fractional derivatives. In particular, we construct polynomial un 
such that 


RL x V{u N {x k ) = RL x V{u{x k ) = f(x k ), ujv(1) = «(1) = 0, 

k = 1,2, • ■ ■ ,N; — 1 < xi < ■ ■ - < Xn < 1, (4.1) 

which is equivalent to find the solution for the following factional initial value problem 
of order s £ (0,1), 

RL x V[u(x) = f(x), u(l)=0, x £ [—1,1]. (4.2) 

4.1. Petrov-Galerkin spectral method for (14.2p . Following the same method¬ 
ology used in [3] for solving (14.21) . to find My £ + J 7< ^ S ' ^ such that mat( 1) = 0 and 

( RL x W 1 u N ,v) = (f,v),VveV N . 


(4.3) 
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Fig. 3.2. Curves of ^(u — un) for Example 3.1 using Gauss-Lobatto points, where 
denotes the corresponding superconvergence points. 



Fig. 3.3. Curves of (u — un) for Example 3.1 using left Gauss-Radau points, where 

denotes the corresponding superconvergence points. 















































Super convergence of fractional polynomial interpolation 


9 


We write the numerical solution as an expansion of GJF basis 

N 

u N (x) = J2u { n s)+ Ji~ s ’- s \x) e + 4r ,-5) (A), (4.4) 

n =0 

Taking vj v = in Eq. (14.3D and using the orthogonality of Legendre polyno¬ 

mials, we derive from m that 

/ \ 77,1 ~ 

“» = r^TTTT T)‘- 0 5 s N - (4 ' 5) 

where {/„} is the sequence of coefficients for the Legendre expansion of function /( x). 
Substituting (14.51) into (14.41) . we obtain the numerical solution: 

N i 

U N (x) = ^ . * ' /rt + (x). (4.6) 

^ o r(n + s + i) 

We would like to indicate the relation of the spectral collocation method (4.1) 
and the Petrov-Galerkin method (4.3). If we select collocation points Xk in (4.1) 
as the Gauss points, i.e., the roots of the Legendre polynomial of degree N, and 
construct an interpolant In f £ Vn such that (Inu)( 1) = 0 and ( lNf){xk ) = f{xk), 
for k = 1, 2,..., N. Denote Wk as weights of the V-point Gauss quadrature. Next, 
we multiply both sides of (4.1) by Lj(xk)wk, sum up over k , and obtain 

(■ RL x V\u n , Lj) = (. I N f , Lj),Vj = 1,2,..., N - 1, (4.7) 

( RL x V>n, L N y = (. I N f , L N )*. (4.8) 

Here * indicates that the integration is carry out by the numerical quadrature. Note 
that the IV-point Gauss quadrature is exact for polynomials of degree 2 N — 1. We see 
that (4.1) is “almost” equivalent to (4.3) in that only one term (14.81) is done by the 
IV-point Gauss quadrature. 

Now, we recall the error estimate of fractional derivative in L 2 -norm as follows. 
Theorem 4.1. (see J3$)Let u and un be the solution of (14.21) and (14.31) . respec¬ 
tively. If f £ C(A) and f^ £ L 2 (A) for all 1 < l < m, then we have that for 

1 < m < N+ 1, 


II RL X V[(U - Ujv)|| + ||u - U N \\ < cA/'- T "||/^ 7 "- ) || w (m-l,m-l) , 

where c is a positive constant independent of u, N and m. 

The above error estimate of fractional derivative leads to also the approximation 
error of function / using Legendre polynomials. Thus, it naturally yields to the 
following superconvergence points for fractional derivatives. 

Proposition 4.2. For the Petrov-Galerkin spectral method for (14.21) . the left 
fractional Riemann-Liouville derivative of numerical solution R1 yDluN superconverges 
to RL x V\u at the zeros of Legendre polynomial P^fl{x), namely, the points {C^} 
satisfying 


Pn+ i(Cf 0) ) = 0, k = 1, 2, - ■ ■ , N + 1, 
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which are usually called Gauss points. 

Due to the special design of basis functions based on different fractional deriva¬ 
tives, the resulting superconvergence points are the same Gauss points for all cases. 

From (14.51) we see that the resulting system is diagonal, which indicates the direct 
expansion of the numerical solution using | + Jr, s ’ s ^(x)|. Therefore, the supercon¬ 
vergence phenomena can also be observed as follows. 

Proposition 4.3. For the Petrov-Galerkin spectral method for m , the nu¬ 
merical solution un superconverges to the exact solution at the zeros of + J^+[ s (a;), 
namely, the points { 77 ^} satisfying 

+ J ( N+r\rii s) ) = 0, k = 1 , 2 , ■ ■ ■ ,N + 1 . 

In the above proposition, superconvergence points { 77 ^}, for the function value, 
depends on the order of fractional derivative. On the contrary, Proposition ^. 21 tells us 
that the fractional derivative of the approximation superconverges at Gauss points, 
which are s- independent, for any 0 < s < 1. The properties of the two kinds of 
superconvergence points manifest the essential differences of the estimates in Theorem 

ro 

The following remark gives us the opportunity to see further how the GJFs of 
(12.51) and (12. 6JI becoming excellent candidates for solving the fractional equation with 
underlying singularities at the boundaries. 

Remark 4.4. For solving problem (14.21) . GJF-Petrov-Galerkin spectral method 
with m as basis functions are applied. While, for the following initial value problem 
involving left Riemann-Liouville fractional derivative: 

RL 1 V s x u(x) = f(x), u(-l) = 0, a: €[-1,1]. (4.9) 

one use m as basis to construct Petrov-Galerkin spectral method, and the main 
results of error estimates and superconvergence points could also be proved. 

In the last part of this section, we present numerical examples to confirm the 
theoretical results and show the superconvergence points for the function value and 
fractional derivatives. 

4.2. Numerical Examples. In this subsection, we report on numerical results 
for problem (14.21) with three different cases of /. 

Example 4.1 Take f(x) = 1 + x + cos(x ) + sin(a;) in (14.21) . 

Example 4.2 Take f{x) = e sln ( x )+ 2 in ()4.2D . 

Example 4.3 Take f(x) = (1 + a:) 7 ' 89 in (14.21) . 

From Theorem 14.11 we know that the errors of function value and its fractional 
derivatives decay exponentially regardless of the unknown solution is singular at 
boundaries or not. Here, we focus on displaying the errors and the superconvergence 
points for these two cases. 

We compute a reference exact solution with N = 41 in (14.61) . Maximum errors 
between the exact solution and the numerical solution with N = 9 for the above 
examples are demonstrated in Figures 14.1114.31 For different orders of fractional 
derivatives as s = 0.1,0.3,0.55,0.7,0.9, the superconvergence points are marked in 
different colors of asterisks respectively. It is easily seen that the errors at these points 
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are much smaller than the global maximum errors. An observation worth noting is 
that when the order s —> 1, the amplitudes of the error curves getting smaller, each 
series of the superconvergence points seems converge to the standard case of s = 1. 

For the purpose of comparison, we present fractional derivative errors in Figures 
14.4114.61 for three examples. In contrast to the function value approximation, super¬ 
convergence points are the same for different orders of fractional derivatives, which 
confirms the conclusion drawn in Proposition 14.21 



Fig. 4.1. Curves of u — with N=9 for Example f.l. 

Remark 4.5. We would like to indicate that although our analysis is for a very 
special case (4.2), our findings of superconvergence points are actually valid for more 
general situations as evidence by the following example. Consider 

R X D 1 u(x)+u(x) = f(x), u( 1) = 0, x e [-1,1]. (4.10) 

We apply the Petrov-Galerkin spectral method (as described above) with N = 9 for 
the exact solution u(x) = (1 — x) 12+a . The error curves of function values for different 
orders of fractional derivative are displayed in Figure 14771 where stars mark the same 
points as in Proposition 14.31 We notice that superconvergence points are the same as 
for problem (14.1011 . 

5. Concluding remarks. In this paper, we discuss superconvergence phenom¬ 
ena for two kinds of spectral interpolations involving fractional derivatives. The in¬ 
tended application is the development of high-order methods for fractional problems 
using spectral methods. 

When interpolating function values, superconvergence points for the Legendre 
basis are located by error analysis. This analysis unifies the identification of super¬ 
convergence points for fractional and first-order derivatives. When interpolating the 
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Fig. 4.2. Curves of u — ujy with N=9 for Example f.2. 



Fig. 4.3. Curves of u — ujy with N=9 for Example f.3. 
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Fig. 4.4. Curves of RR T>f(u — un) with N=12 for Example f.l. 



Fig. 4.5. Curves of RR X>f(u — ujv ) with N=18 for Example f.2. 
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Fig. 4.6. Curves of RR T>f(u — un) with N=10 for Example f.3. 



Fig. 4.7. Curves of u — un with N=9 for Eq. (14. 1 01) 
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fractional derivative values, we found significant differences between super convergence 
points for function values and fractional derivatives. Numerical comparisons between 
the two cases are provided. 

In the future, we plan to study more general fractional differential equations such 
as ([4.10ft . 
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